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Abstract 

Under panicky situation, human have tendency to rush toward a 
particular direction for escape. I show here that this tendency alone 
causes increase in crowding and which could eventually trigger jam- 
ming that is not preferable. Further, it is proposed that potential flow 
^^ ■ theory can be employed in finding suitable strategy for escape. 

No one would prefer to be present in unstable situation of over-crowding 
leading to jamming and stampedes when individuals try to escape the sit- 
uation under panic. Unfortunately, panicky situations are not completely 
avoidable as they are caused by forces, e.g. fires, whose origins in space and 
time remain unpredictable. We could only hope that spaces used by human 
and their escape routes were well laid out to avoid any casualty if panicky 
situation arises. The designs of layouts and appropriate escape routes require 
knowledge and mathematical representation of individual and collective be- 
haviors of human during the escape. Based on the knowledge, we could 
adopt mainly two kinds of design strategies. In the first strategy (referred 
to as SI hereafter), spaces and escape routes can be designed compatible 
to the natural behavior of human and without implementing any ways to 
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change or manipulate human behavior during escape. This kind of strategy 
is usually employed in practice (see Ref. [1] and references cited therein). In 
the second strategy (referred to as 52 hereafter), the design can incorporate 
certain ways to change human behavior and their motions during the escape. 
Here, first I show that the natural behavior of human would always lead 
to over-crowding in 51 strategy and then I describe S2 strategy suggesting 
potential flow theory based motions for individuals to avoid or minimize the 
over-crowding. 

I put myself, in a hypothetical situation, among the people which were 
dancing and were dispersed uniformly in a large hall. Suddenly music stopped, 
there was a panic, reason unknown, and individuals started rushing directly 
toward the entrance, also an only exit for the hall. In the beginning of panic, 
I was able to move freely toward the exit with velocity v^(t)e^(t). But soon I 
realized that I had to be careful not to bump into someone and my behavior 
was implicitly controlled by surrounding individuals and crowd. I started to 
have difficulty moving directly toward the exit and continued moving with 
velocity v^(t) where ever I could, but was always trying hard to reach the 
exit. I was not fortunate enough and soon afterward got caught in the jam 
situation when the crowd was almost stand still and pushing, each other des- 
perately, to squeeze through their way to the exit. People did not remain 
gentle then and what happened afterward is left for imagination. 

Though I was exhausted, I did realize that I experienced and felt the 
mathematical equation suggested by Helbing et al. [1] 

ro> d^) = m <ffle?W-v, W + ^ ^ + ^ ^ (1) 

Here mi is ith individual's mass (e.g. my mass), v®{t)e®{t) is intended ve- 
locity of an individual toward the exit due to natural tendency of human, 
Vj(t) is actual velocity of an individual with a certain characteristic time 
Tj and J2j(=£i) fij is sum of anticipated forces to keep myself a safe distance 
(if possible) from other surrounding individuals j's and contact forces dur- 
ing pushing. Also, J2 W ftw is anticipated force to keep myself away from the 
nearest walls, had I been caught between the crowd and walls of the hall. 
The Lagrangian Eq. (1) is the most reasonable mathematical description, 
available so far, for individual's behavior under the influence of surrounding 
individuals' behavior and presence of walls in panicky situation. Helbing et 



al. [1] considered a few benchmarks to show the robustness and usefulness 
of this equation for designing escape strategy of type SI. 

So for strategy SI, now I show that the behavior governed by Eq. (1) 
would always result in over-crowding if every individual in the hall rushes 
directly toward the exit. For the discussion, consider magnitude v®(t) of the 
intended velocity to be identical for all individuals, independent of time and 
equal to v°. But the unit direction vector e°(t) would be different depending 
on the position of the individual with respect to the exit. Also, consider floor 
of the hall in x — y plane of a rectangular coordinate system x — y and exit 
of the hall at the origin. Every individual's position can then be given by 
coordinates (x, y) and thus the intended velocity, directly toward the origin, 
of any individual located at (x, y) can be given by the Eulerian velocity field 
v°(x, y,t), written as 

v°(x,y,t) = v°[- r^—4 ~ T^T^ j] ( 2 ) 

L v x 2 + y 2 ^/x 2 + y 2 J 

where i and j are unit vector along the x and y axes, respectively. Using 
the usual recipe of Kinetic Theory approach, Eulerian equations for number 
density N(x, y, t) of people and number weighted velocity field ~V(x, y, t) can 
be obtained by using Eq. (1) and Lagrangian equation — -^- = Vj(i) for 
position vector R«(t) of each individual. The governing equation for N and 
V can be written as, respectively, 

R N = [JL + V • V] JV = -NdivV, (3) 

V(x, y, t) = w\x, y, t) - r^V + V c , (4) 

where v°(x, y,t) is any intended velocity field (e.g. as given by Eq. 2) and 
V c accounts for the velocity due to forces f^ and f iw . While writing Eq. (4), 
time scale Tj for all individuals is assumed to be identical. The three terms 
on the right hand side (rhs) of Eq. (4) can be considered as zeroth order, 
first order and second order velocity terms, respectively, based on the order 
in which corresponding events causing these different velocities occur. 

In the beginning of the panicky situation (Stage 1) when I was able to 
move freely, V c = and I immediately gained intended velocity, i.e. 

V(x,y,t)^v°(x,y,t), (5) 



as suggested by Eq. (4). In the beginning of Stage II, when I was avoiding 
collision with other individuals and was adjusting myself with small response 
or characteristic time scale Tj, V c was negligible. For this stage, Eq. (4) 
along with perturbation expansion around v°(x,y,t) suggest 

V(x, y, t) = v°(x, y, t) - r^ + v° • V]v°. (6) 

These velocity fields of initial two stages are sufficient to show increase in 
over-crowding or in number density. 

Now during Stage I, Eqs. (2) and (5) provide divY = —v°/y/x 2 + y 2 and 
its substitution into Eq. (3) yields 



DN/Dt = v°N/^x 2 + y 2 >0, (7) 

suggesting increase in number density iV in time. This increase continues 
in time during the beginning of the Stage II It should be noted that an 
additional Ti containing term on the rhs of Eq. (6) can not reduce or increase 
the rate of increase, as div [(J^ + v° • V)v°] = when Eq. (2) is used for v°. 
This exhibits clearly that individuals' intentions to rush directly toward the 
exit is a main cause for increase in time of the number density of people i.e. 
over-crowding. 

Now I describe strategy 5*2. The cause of over-crowding is due to the com- 
pressible nature of the velocity field V, i.e. non-zero value for divV . The 
intended velocity field v° arising from human's natural tendency is compress- 
ible too. So in order to avoid over-crowding and subsequent triggered events, 
the first required condition is that velocity field V should be incompressible 
i.e. divV = 0. Further, during escape I would not prefer to be moving in a 
circle again and again inside the hall. Thus, the second required condition 
of zero circulation suggests curlV = 0. These two conditions suggest the 
velocity field to be identical to that of inviscid, incompressible, irrotational 
fluid flow or potential flow, written in terms of a velocity potential function 
(j) as 

V = grad</), V 2 = 0. (8) 

Thus, if individuals move along the streamlines of appropriate potential flow 
with suitable speed, triggering of over-crowding and subsequent disastrous 
events can be altogether avoided. There is a possibility that a few individuals 
may not move exactly as required and perturbation in their behaviors arise 
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due to finite response time to follow streamlines. The strategy £2 can then 
be further refined by using strategy SI. To do that, potential velocity field 
can be considered as intended velocity v° in the framework of SI. The over- 
crowding will be caused at first order by the compressible nature of the second 
term on the rhs of Eq. (6). So for potential velocity field v°, divergence of 
Eq. (6) yields 

< 0, (9) 

where each term on the rhs is proportional to square of different strain rate. 
The Eqs. (9) and (3) suggest that over-crowding is likely to occur in the 
region of high strain rate of potential velocity field. Though it is not possible 
to have these strain rates zero everywhere, spaces and escape routes should be 
designed in a manner to minimize, as far as possible, areas of high strain rates 
in the potential velocity field. In addition, the quick response of human, i.e. 
small value for Tj, would be an added advantage to reduce the over-crowding. 
By now, we know suitable strategy. Now, important questions arise: 
How to implement S2 strategy in reality? How to manipulate individuals' 
behaviors for their own advantage during escape so they quickly respond 
and move in, already established, streamlines pattern in panicky situations? 
I leave these difficult tasks for my fellow engineers to come up with clever 
ways, including and other than, marking streamlines of appropriate potential 
flow on the floor, using moving/flashing lights for speed at every location to 
guide individuals during escape. 
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